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^ , The anomalously large cos(2</)) asymmetry measured in the Drell-Yan process is 

discussed. Possible origins of this large deviation from the Lam- Tung relation are 
considered with emphasis on the comparison of two particular proposals: one that 
' suggests it arises from a QCD vacuum effect and one that suggests it is a hadronic 

effect. Experimental signatures distinguishing these effects are discussed. 
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1. Introduction 



Qh! Azimuthal asymmetries in the unpolarized Drell-Yan (DY) process differ- 

^ • ential cross section arise only in the following way 

(U ! ^ ^ 

■ ^ oc (l + Acos^e + /isin26' cos(/) -f -sin^ 6'cos20) , (1) 

cr ail V 2 / 

where is the angle between the lepton and hadron planes in the lepton 

rS \ center of mass frame (see Fig. 3 of Ref.-*^). In the parton model (order a^) 

' quark-antiquark annihilation yields \ = \, ^ = v = Q. The leading order 

(LO) perturbative QCD corrections (order a].) lead to ^ 0, v ^ Q and 

A 1, such that the so-called Lam- Tung relation 1 — A — 2i^ = holds. 

Beyond LO, small deviations from the Lam- Tung relation will arise. If one 

defines the quantity k = —\{^ — A — 2v) as a measure of the deviation 

from the Lam- Tung relation, it has been calculatea^ ^ that at order a1 k 

is small and negative: — k^O.01, for values of the muon pair's transverse 

momentum Qt of up to 3 GeV/c. 

Surprisingly, the data is incompatible with the Lam- Tung relation and 

with its small order-a? modification as welP. These data from CERN's 

NAIO CoUaboratioi^EIand Fermilab's E615 CoUaboratiori^ are for tt'TV ^ 

fi^fi~X, with N ^ D and W. The 7r~-beam energies range from 140 GeV 
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up to 286 GeV and the invariant mass Q of the lepton pair is in the range 
Q 4—12 GeV. The measured values for k are an order of magnitude 
larger than the order-a^ result and moreover, of opposite sign. 

Several explanations have been put forward, but not all of them will 
be reviewed here. Some unlikely explanations would be: i) NNLO pQCD 
corrections could solve the discrepancy (but in that case the perturbative 
expansion itself would be questionable); ii) it could be a higher twist effect 
(but Q"^ > 16 GeV^ seems too high and according to the Fermilab data the 
deviation disappears at high x^, contrary to higher twist expectation; also, 
one would expect fi > v, whereas in the data v ^ fj, ~ 0); in) it could be a 
nuclear effect, since (j(Qt)w / <^{.Qt)d is an increasing function of Qt (but 
according to Ylei^ v{Qt) shows no apparent nuclear dependence). 

The two possible explanations that will be discussed and compared here 
are: i) a QCD vacuum effec^ ii) a hadronic effect, arising from non- 
collinear parton configurations^ The following will largely be based on a 
recent comparative study performed in collaboration with A. Brandenburg, 
O. Nachtmann and A. UtermanrH 

2. Explanation in terms of a QCD vacuum effect 

Usually the DY process at Q ^ 4 — 12 GeV is described by collinear fac- 
torization. CoUinear quarks inside unpolarized hadrons are unpolarized 
themselves, implying a trivial quark-antiquark spin density matrix: 



The QCD vacuum may alter this. The gluon condensate leads to a chro- 
momagnetic field strength (Savvidy; Shifman, Vainshtein, Zakharov; ...) 



with giuon fields having a typical correlation length a k, 0.35 fm in Eu- 
clidean space. Taking this to be an invariant length in Minkowski spacJ^ 
leads to the picture of a fluctuating domain structure of the vacuum with 
typical domain size a, schematically depicted in Fig. ^ If a fast hadron, and 
with it a fast quark, traverses this domain structure, the time for traversing 
a vacuum domain is of the order of the correlation length: t « a. Due to 
the presence of a background chromomagnetic field the quark will acquire 
a transverse polarization (the Soko loy-T ernov effect). The time to build up 
transverse polarization is estimatedSEl to be much shorter than the time it 
takes to traverse the domain, i.e. t <^ a. The radiated gluons/photons are 



P' 




(2) 



(g'^B'^ix) ■ B^ix)) « (700MeV)^ 



(3) 
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Figure 1. Left figure: cartoon of the chromomagnetic field domain structure of the 
QCD vacuum. Right figure: a fast quark traversing a domain. 

just part of the cloud of virtual particles; in other words, they are included 
in the wave function. There will be no average polarization. However, if the 
quark will annihilate with an antiquark in a high energy scattering experi- 
ment, such as DY, the polarization of the quark and the antiquark may be 
correlated if they annihilate within a certain domain. Therefore, the QCD 




Figure 2. Annihilation of a quark and an antiquark inside the same domain. 

vacuum can induce a spin correlation between an annihilating q q pair. The 
quark-antiquark spin density matrix Eq. jSJ will then be modified into 

p^'''^'^ = + FjeTj (»1 + Gjl<» CTj + H^j (Ti ® CTj} . (4) 

Only if Hij — FiGj, then the spin density matrix factorizcs. But this is not 
necessarily so, in which case it could be called entangled. Brandenburg, 
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Nachtmann & MirkeJ^ demonstrated that the diagonal elements Hn and 
H22 can give rise to a deviation from the Lam- Tung relation: 

A simple assumption for the transverse momentum dependence of {H22 ~ 
Hii)/{1 + H33) produced a good fit to the data: 

= '^o — — — 4- , with Ko = 0.17 and ttit = 1.5 GeV. (6) 

Note that for this Ansatz k approaches a constant value (kq) for large Qt- 
In other words, the vacuum effect could persist out to large values of Qt- 
The dependence of the vacuum effect is not known, but there is also no 
reason to assume that the spin correlation due to the QCD vacuum effect 
has to decrease with increasing . 



3. Explanation as a hadronic effect 

Usually if one assumes that factorization of soft and hard energy scales in 
a hard scattering process occurs, one implicitly also assumes factorization 
of the spin density matrix. In the present section this will indeed be as- 
sumed, but another common assumption will be dropped, namely that of 
collinear factorization. It will be investigated what happens if one allows for 
transverse momentum dependent parton distributions (TMDs). The spin 
density matrix of a noncoUinear quark inside an unpolarized hadron can 
be nontrivial. In other words, the transverse polarization of a noncoUinear 
quark inside an unpolarized hadron in principle can ha ve a preferred direc- 
tion and the TMD describing that situation is called hi li^. As pointed out 
in Rcf !^ nonzero hi leads to a deviation from Lam- Tung relation. It offers 
a parton model explanation of the DY data (i.e. with A = 1 and ^ = 0): 
K — ^ oc hi (tt) hi (N) . In this way a good fit to data was obtained 
by assuming Gaussian transverse momentum dependence. The reason for 
this choice of transverse momentum dependence is that in order to be con- 
sistent with the factorization of the cross section in terms of TMDs, the 
transverse momentum of partons should not introduce another large scale. 
Therefore, explaining the Lam- Tung relation within this framework neces- 
sarily implies that k = ^ for large Qt- This offers a possible way to 
distinguish between the hadronic effect and the QCD vacuum effect. 

It may be good to mention that not only a fit of hj^ to data has been 
made (under certain assumptions), also several model calculations of h^ 
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and some of its resulting asymmetries have been performec^^^^l, based 
on the recent insight that T-odd TMDs hke arise from the gauge hnk. 

In order to see the parton model expectation k = ^ ^ at large Qt in 
the data, one has to keep in mind that the pQCD contributions (that grow 
as Qt increases) will have to be subtracted. For k perturbative corrections 
arise at order aj., but for v already at order as- To be spe cific, at large Qt 
hard gluon radiation (to first order in Us) gives rise tcM^ 

Due to this growing large-Qr perturbative contribution the fall-off of the 
hi contribution will not be visible directly from the behavior of ly at large 
Qt- Therefore, in order to use v as function of Qt to differentiate between 
effects, it is necessary to subtract the calculable pQCD contributions. In 
Fig. 01 an illustration of this point is given. The dashed curve corresponds 




Figure 3. Impression of possible contributions to u as function of Qt compared to DY 
data of NAIO (for Q = 8 GeV). Dashed curve: contribution from perturbative one-gluon 
radiation. Dotted curve: contribution from a nonzero h^. Solid curve: their sum. 



to the contribution of Eq. ^ a,t Q — 8 GcV. The dotted line is a pos- 
sible, parton model level, contribution from with Gaussian transverse 
momentum dependence. Together these contributions yield the solid curve 
(although strictly speaking it is not the case that one can simply add them, 
since one is a noncollinear parton model contribution expected to be valid 
for small Qt and the other is an order-cts result within coUinear factor- 
ization expected to be valid at large Qt)- The data are from the NAIO 
Collaboration for a pion beam energy of 194 GeV/cl^. 

The Q^ dependence of the contribution is not known to date. Only 
the effect of resummation of soft gluon radiation on the contribution to 
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V (a nd K ) has been studied to some extent and was found to be quite impor- 
tant^!. It gives rise to a considerable Sudakov suppression with increasing 
Q: in going from Q = 10 to 90 GeV, the contribution decreases by an order 
of magnitude and approximately follows a 1/Q behavior (although it is nei- 
ther a dynamical nor a kinematical higher twist effect). Interestingly, the 
contribution from hard gluon radiation Q decreases more rapidly: as l/Q^ 
at fixed Qx- But it seems safe to conclude that using the dependence 
of V (or k) to differentiate between effects is not feasible at present. 

By assumption, nonzero h-^ gives rise to a factorized product of spin 
density matrices p^"?'?) = p('J) (55 p^?) witf!^ 

Therefore, 7?^ = FiG^ with TJas = 0. Unfortunately it is hard to observe 
the difference between iJaa = and ^^33 ^ 0. But the factorization Hij = 
FiGj should shows itself via consistency among various processes, which is 
based on the fact that the same function hi appears in different proces ses. 
Regarding this universality, complications have recently been addressec^^ 
that go beyond the sign change that occurs between semi-inclusive DIS 
[ep e' TT X) and DY: (/i5'")sidis — ~{hi)DY- Nevertheless, the different 
numerical factors with which hi arises in different processes are calculable 
(functions of Nc only) and can be taken into account. 

4. Hadronic effect versus vacuum effect 

Summarizing the features of the two approaches in a table: 



Tabic 1. Comparison of the hadronic and the vacuum effect 





hi 


7^0 


QCD vacuum effect 




pM 




possibly entangled 


Q dependence 


K ~ 


1/Q 


? 


Qt 00 


K — 


» 


need not disappear (k — > kq) 


flavor dependence 


yes 




flavor blind 


X dependence 


yes 




yes, but flavor blind 



As indicated in the table, the hadronic effect will generally be flavor depen- 
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dent and have an x dependence that is flavor dependent, since there is no 
reason to assume that hj^ for the u quark should be the same as (or simply 
related to) that for the d quark. This is different from the QCD vacuum 
effect, which in this sense is flavor blind; it does not matter whether the 
spin correlation is between uu or dd (except for presumably small mass 
corrections). There will be an x dependence, since that determines the 
energy of the annihilation process, but this again should be flavor blind. 
It should be emphasized that flavor blindness in general does not imply 
hadron blindness or even process blindness. So the best next step would be 
to perform experiments with different beams (tt^ ,p,p, . . where 7r+ and 
p offer the advantage of having valence anti-quarks) and in different kine- 
matical regimes. For instance, the measurement of (cos 20) can be done at 
RHIC in pp ii^fi^X, or in pp — » ^~^ii~X at Fermilab or GSI/FAIR. 

The use of polarized beams can also help (e.g. at RHIC or GSI). In the 
DY process with one transversely polarized hadron, the differential cross 
section can namely depend on the azimuthal angle (f>s of the transverse 
hadron spin (St) compared to the lepton plane: 



dQ d(f>s 



cx 1 + cos^ 9 + sin'^ ( 



1/ 



— cos 20 ~ p \Srp\ sin(0 + ^s) 



Within the framework of TMDs the analyzing power p is proportional to 
the product hj^ hi ^, which involves the transversity function hi. A nonzero 
function hi will provide a relation between v and p, which in case of one 
(dominant) flavor (usually called w-quark dominance) and Gaussian trans- 
verse momentum dependences, is approximately given by 



^ r^^' (10) 



2 V i^max /l ' 

where I'max is the maximum value attained by v{Qt)- This r elation de- 
pends on the magnitude of hi compared to /i (see Refs.^^l for explicit 
examples) and this may be extracted from double transverse spin asym- 
metries in DY (potentially at RHIC or GSI) or from SIDIS data (from e.g. 
HERMES or COMPASS) by exploiting the interference fragmentation func- 
tions (which can be obtained from e+e" data, e.g. at BELLE). 

Also semi-inclusive DIS can be used. The (cos2(/)) in ep e' tt X 
would be cx hiH^, where is the Collins fragmentation function (also 
obtainable from BELLE) . This particular SIDIS observable has been studied 
using model calculation^^. All this illustrates how the consistency among 
processes may be used to test the hi hypothesis. 
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5. Conclusions 

A transverse spin correlation in quark- antiquark annihilation {q^q^ 7*) 
will lead to a cos(2(/)) asymmetry in the DY lepton-pair angular distribution. 
Such a spin correlation can arise from the chromomagnetic background field 
in the QCD vacuum or from noncoUinear partons. If a flavor dependence is 
observed in future data, it would favor a hadronic effect. On the other hand, 
persistence of the asymmetry at large values of Qt and Q (after subtraction 
of pQCD corrections if needed) would favor a vacuum effect. Several future 
and ongoing experiments will be able to provide crucial information on 
these dependences. 
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